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ABSTRACT 

The integrated ultraviolet flux arising from QSOs and/or hot, massive stars in metal-producing 
young galaxies is likely responsible for maintaining the intergalactic diffuse gas and the Lya forest 
clouds in a highly ionized state. The spectrum and intensity of such UV background have generally 
been obtained by modeling the reprocessing due to intervening material as a pure photoelectric 
absorption process. However, owing to the emission from radiative recombinations within the ab- 
sorbing clouds, a photoionized clumpy medium could contribute substantially to the metagalactic 
flux. In other words, QSO absorption-line systems are sources not just sinks of ionizing photons. 

We present a detailed calculation of the propagation of AGN-like ionizing radiation through 
the intergalactic space. We model the ionization state of absorbing clouds, and show that the 
universe will be more opaque above 4 Ryd than previously estimated. Singly ionized helium 
in Lya forest clouds and Lyman-limit systems is found to be very efficient in reprocessing soft 
X-ray, helium-ionizing photons into ultraviolet, hydrogen-ionizing ones. We demonstrate that 
a signiflcant fraction of the absorbed primary photons (emitted, e.g., by quasar sources) will 
be reradiated by the photoionized gas through Lya line emission, two-photon continuum, and 
recombination continuum radiation. 

In the light of new data and recent studies, we also reassess the contribution of the QSOs 
observed in optical surveys to the UV extragalactic background, and flnd that the stochastic 
reprocessing of quasar Lyman continuum radiation by hydrogen and helium along the line of 
sight will signiflcantly affect the amplitude, spectral shape, and fluctuations properties of the 
metagalactic flux. In a scenario in which QSOs are the primary source of ionizing photons in the 
universe, the integrated H I Lya emission at 2; = from photoionized Lya clouds and Lyman- 
limit systems is found to be at a level of less than 5% of current observational limits on the far-UV 
extragalactic radiation flux. We show that J912 increases from 10~^^ ergscm"^ Hz""^ sr~"^ 
at the present epoch to 5 X 10~^^ ergscm"^ Hz""^ sr~"^ at z = 2.5. The attenuated direct 
quasar emission plus recombination radiation from intergalactic gas appears to provide enough 
hydrogen-ionizing photons to satisfy the proximity effect at large redshift. The He II /H I ratio 
in the diffuse intergalactic medium and the Lya clouds increases from 25 at 2; = to 45 at 
z = 2.5, to decrease again below 30 for z > 4.5. The spectrum of the ionizing background at high 
redshift is shown to have a hump at energies below 40.8 eV due to redshift-smeared He II Lya line 
and two-photon continuum emission. We propose that observations of low-ionization species such 



as O II in metal-line absorption systems may be able to test the presence of such a prominent 
feature in the UV background spectrum. We also note that, if the metagalactic flux is dominated 
by QSOs, the suggested steep decline of their ionizing emissivity beyond 2; ~ 4 should produce 
an increase in the observed rate of incidence of Lya forest clouds at these redshifts relative to an 
extrapolation from the intermediate- 2; regime, as observed by Williger et al. (1994). 

Subject headings: cosmology: observations - diffuse radiation - intergalactic medium - quasars: 
absorption lines 







1. INTRODUCTION 

The existence of a uniformly distributed intergalactic medium (IGM), which contains the 
bulk of the baryonic matter of the universe in most models of structure formation, is predicted as 
a product of primordial nucleosynthesis. The hydrogen component of this IGM must have been 
highly ionized by 2; ~ 5, as the flux decrement observed on the blue side of the Lya emission 
line of neutral hydrogen in the spectra of high redshift QSOs appears entirely consistent with the 
blanketing by discrete absorption lines along the line of sight (Steidel & Sargent 1987; Giallongo 
et al. 1994). It is widely believed (but see, e.g., Sciama 1993) that the integrated ultraviolet flux 
arising from QSOs and/or hot, massive stars in metal-producing young galaxies is responsible 
for maintaining the intergalactic diffuse gas and the Lya forest clouds in a highly ionized state. 
Such UV background may also be responsible for the ionization of the metal-rich QSO absorption 
systems (Steidel 1990; Vogel & Reimers 1993) and of the hydrogen clouds located in the galactic 
halo (Fransson & Chevalier 1985; Ferrara & Field 1994), and for producing the sharp edges of 
H I disks observed in nearby spiral galaxies (Bochkarev & Sunyaev 1977; Maloney 1993; Corbelli 
& Salpeter 1993; Dove & ShuU 1994). Photoionization by UV radiation may also inhibit the 
collapse of small-mass galaxies (Dekel & Rees 1987; Efstathiou 1992). 

The Hubble Space Telescope Faint Object Camera (FOG) detection of redshifted He II 304 A 
absorption in the spectrum of the z = 3.29 quasar Q0302-003 (Jakobsen et al. 1994) has recently 
provided new constraints on the thermal and ionization state of the IGM, and on the spectrum 
of the UV metagalactic flux at early epochs. The high ratio of He II to H I which is necessary to 
account for the observed absorption trough requires a soft radiation fleld at 4 Ryd compared to 1 
Ryd (Madau & Meiksin 1994; Giroux, Fardal, & ShuU 1995). Although the mean spectrum of the 
individual sources which dominate the metagalactic flux might not be very steep, the cosmological 
"flltering" through material along the line of sight is known to signiflcantly soften the background 
intensity. Absorption by intervening matter actually makes the universe at high redshift optically 
thick to Lyman continuum (LyG) photons: enough neutral hydrogen and singly ionized helium is 
contained in the intergalactic clumps of highly ionized gas which form the Lya forest, and in the 
metal-line absorption systems associated with the halo regions of bright galaxies, to signiflcantly 
attenuate the LyG flux from all but the nearby sources (Bechtold et al. 1987; Miralda-Escude 
& Ostriker 1990; M0ller & Jakobsen 1990; Madau 1991, 1992; Meiksin & Madau 1993; Zuo & 
Phinney 1993; Madau 1995). 

This is the second paper in a series aimed at a detailed numerical study of the absorption 
and reprocessing of UV photons in a clumpy, photoionized universe. In Paper I (Madau 1995), 
we have assessed the effects of the stochastic attenuation produced by QSO absorption systems 
along the line of sight on the broadband colors of galaxies at cosmological distances. It is the 
purpose of this paper to focus on the effects of intervening discrete absorbers on the propagation 
of ionizing diffuse radiation. The primary motivation of this study is the realization that the 
integrated UV background will actually be the sum of the direct radiation from discrete sources 
and the emission from radiative recombinations within the absorbing clouds and the intercloud, 
diffuse IGM. Since the intercloud medium is optically thin to LyG photons, its contribution to the 
photoionizing radiation fleld will be negligible. We show in this paper that this will not be true, 
however, for the cloudy component, as a signiflcant fraction of the absorbed primary photons will 
be reradiated by the photoionized gas through Lya line emission, two-photon continuum, and 
recombination continuum radiation. The plan is as follows. In § 2 we review the basic theory 
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of cosmological radiative transfer in a clumpy universe. In § 3 we model the physical state of 
absorbing clouds, as the knowledge of their ionization structure is critical to understanding the 
propagation of ionizing photons through intergalactic space. In § 4 we determine the contribution 
of recombination radiation from Lya forest clouds and Lyman-limit systems to the metagalactic 
flux. In § 5 we numerically compute the spectrum of the background radiation fleld as a function 
of redshift, assuming that the quasars detected in optical surveys are the primary source of LyC 
photons in the universe. We derive expressions for the hydrogen and helium ionization rates, 
and discuss the implications of our results for the Gunn-Peterson test, the proximity effect, the 
He II Lya cosmic opacity, the fluctuations of the metagalactic ionizing flux at high redshift, and 
the properties of the Lya forest clouds and metal-line absorption systems. We briefly summarize 
our conclusions in § 6. In Appendix A we provide some fltting formulae to the recombination 
rate coefficients, while in Appendix B we present useful analytical approximations to the effective 
cosmic opacity and the recombination emissivity. Finally, in Appendix C we discuss a numerical 
algorithm for a fast iterative solution to the equation of radiative transfer. Readers who are 
interested in the main results and not in the technical details of the model may wish to proceed 
directly to § 5 at this point. 

2. BASIC THEORY 



2.1. Cosmological Radiative Transfer 

The radiative transfer equation in cosmology describes the evolution in time of the specific 
intensity J (in units of ergs cm~^ Hz""^ sr~"^) of a diffuse radiation fleld (e.g., Peebles 1993): 

d a d \ a c 

— - U-— \ J = -3-J - ckJ + —e, (1) 
Ot a ov J a 47r 

where a is the cosmological scale parameter, c the speed of the light, n the continuum absorption 
coefficient per unit length along the line of sight, and e is the proper space- averaged volume 
emissivity. Integrating equation (1) and averaging over all lines of sight yields the mean specific 
intensity of the radiation background at the observed frequency Vg, as seen by an observer at 
redshift Zg'. 



where u = Vo{l + z)l{l + Zo), 



is the line element in a Friedmann cosmology, ifo = 50/i5okms ^ Mpc ^ the Hubble constant, 
and go the deceleration parameter. 

2.2. Intergalactic Absorption 
The effective optical depth T^ff due to discrete absorption systems is defined as 



^eff = -ln((e ^)), 



(4) 
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where (e ^) is the average transmission over aU hnes of sight. For Poisson-distributed clouds we 
have (Paresce, McKee, & Bowyer 1980) 

Teffiiyo,Zo,z) = J dz' dNm ~ g (5) 

where N/ dNjudz' is the redshift and column density distribution of absorbers along the line of 
sight. The continuum optical depth through an individual cloud is (assuming a pure hydrogen- 
helium gas) 

T(iy') = Nm'Jmiiy') + iVHeIfTHel(/^') + iVHeIIfTHeIl(/^'), (6) 

where = i^o(l + -2')/(l + -2o), and cthi , CHel , and (JHell are the hydrogen and helium photoioniza- 
tion cross-sections. When r ^ 1, Tg// becomes equal to the mean optical depth. In the opposite 
limit, the obscuration is picket fence-type, and the effective optical depth becomes equal to the 
mean number of optically thick systems along the line of sight. 

Two classes of absorbers are known to dominate the cosmic LyC opacity, with the Lya forest 
clouds being the most numerous. These have neutral column densities ranging from ~ 10"^^ to 
maybe 10"^^ cm~^, close to primordial metal abundances, and evolve rapidly between 1.7 < z < A 
(Murdoch et al. 1986). The distribution of the number of clouds as a function of H I column 
roughly follows a single power- law dN/dNm oc iVg/, with /3 ~ 1.5 (Tytler et al. 1995a; Tytler 
1987). The redshift evolution of the rarer Lyman-limit systems (log iVni > 17.2), which are 
optically thick to radiation with energy greater than 1 Ryd, is much weaker over the range 
0.7 < z < 3.6 (Storrie-Lombardi et al. 1994; Sargent, Steidel, & Boksenberg 1989). About 2/3 of 
the Lyman-limit absorbers are associated with heavy-element lines. 

To facilitate a comparison with previous estimates of the UV background (Miralda-Escude 
& Ostriker 1990; Madau 1991, 1992; Miralda-Escude & Ostriker 1992; Meiksin & Madau 1993; 
Giroux & Shapiro 1995), we assume in the following a constant value of /3 = 1.5 over 7 decades 
in iVni. For the redshift and column density distribution of intervening absorbers we take 



d^N ^ NoN-l-^(l + zy (10" < iVm < 1.59 X lOi^cm-2); 



dNmdz 1 5.4 x 10^iV£/-^(l + zy-^^ (1.59 x 10^^ < iVm < lO^" cm'^). 



(7) 



where Nq = 2.9 X 10^ and 7 = 2.46 for z > 1.7. With the above normalization, there are 
dN / dz = 4.3(1 -|- zY''^^ Lya forest clouds with rest equivalent width W > 0.32 A (corresponding 
to iVni > 1.8 X 10"^"^ cm~^ for an average Doppler width of 28 km s~"^), as measured at high redshift 
by Press & Rybicki (1993). At 2; ~ 0, the number of Lya forest line is ~ 5 times larger than 
expected from a simple extrapolation of the high redshift evolution (Bahcall et al. 1991; Morris 
et al. 1991). For z < 1.7, we therefore take 7 = 1.0 and normalize imposing continuity at z = 1.7. 
This gives ~ 18 lines per unit redshift with W > 0.32 A at the present epoch, consistent with 
the findings of Bahcall et al. (1993). At high- 2;, the model adopted will provide about the same 
photoelectric opacity as in model A2 of Miralda-Escude & Ostriker (1990), and in the medium 
attenuation model MA of Meiksin & Madau (1993). It is fair to point out that deviation from 
a single power-law in column density could be significant. In particular, there is very limited 
information on Lya forest absorbers with 10"^^ < iVni ^ 10"^^cm~^. These high column density 
clouds will dominate the Lya forest contribution to the effective photoelectric optical depth. 
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making its value rather uncertain. Meiksin & Madau (1993) have suggested that the pubhshed 
data may actuaUy indicate an underdensity of Lya with columns iVni > 10"^^ cm~^. Giallongo et 
al. (1993), Petitjean et al. (1993), and Kulkarni et al. (1995) have also reported evidence for a 
break in the iVni power-law at log iVni ^ 15. This deficit could result in a significant reduction 
of the forest continuum opacity, and therefore in an increased background ionizing flux. We will 
discuss this possibility in more details in § 5. 

In the following, we will also neglect the photoelectric opacity associated with a uniformly 
distributed IGM. This is consistent with the lack of a H I Gunn-Peterson absorption trough in the 
spectra of QSOs at large redshift (e.g., Giallongo et al. 1994). Note, however, that if the recent 
Hubble Space Telescope observation of a He II Lya trough in the spectrum of quasar Q0302-03 
(Jakobsen et al. 1994) were to require a Gunn-Peterson optical depth > 10, the He II LyC opacity 
of a diffuse IGM component would be signiflcant (Madau & Meiksin 1994). 

3. THE IONIZATION STATE OF QUASAR ABSORPTION SYSTEMS 

A self-consistent calculation of the effective optical depth in equation (5) requires a de- 
tailed understanding of the ionization state of intervening material. In particular, the amount of 
He I and He II which is present along the line of sight must be inferred by modeling the phys- 
ical conditions in the absorbing clouds themselves. This will be done in this section. We will 
assume the gas to be in ionization equilibrium with the background UV radiation. Deviation 
from equilibrium will be small, of order the ratio of the ionization timescale to the Hubble time, 
tion/tR ~ 6 X 10~^ Jg^2 -2i[(-'- "I" -2^)/4]^'^^^50, where the background intensity is measured at the 
hydrogen Lyman edge in units of 10"^"^ ergscm"^ Hz~"^ sr~"^. 

In principle, the temperature T of the absorbing gas should also be computed self-consistently, 
by balancing photoionization heating with free-bound and free-free radiative losses, coUisional 
excitation losses, and Compton cooling losses from scattering off the microwave background. The 
temperature fully ionized gas may reach depends, however, on the poorly known cloud density. For 
~ 4 X 10""^ cm~^ , coUisional excitation of neutral hydrogen is the dominant cooling mechanism, 
and the cloud temperature increases with decreasing density, to drop again at lower densities when 
Compton cooling starts to dominate (Meiksin 1994). Moreover, when nn ^ 10""^ cm~^ (the typical 
density of a spherical Lya cloud with neutral hydrogen column 10"^"^ cm~^ and diameter 100 kpc 
in equilibrium with a radiation fleld of intensity J9i2,-2i ~ 1), the cooling timescale is longer than 
the age of the universe at the redshift of interest, and the gas will retain a memory of the initial 
conditions. The contraction or expansion of the gas toward hydrostatic equilibrium will further 
affect its thermal properties. On the other hand, we can estimate the temperature of the absorbing 
gas from the measurements of Lya line proflles in high resolution, high signal-to-noise spectra. 
For simplicity, we shall assume here a flxed temperature ofT = 2.5 X 10"^ K. This is the value 
derived from modeling the broader Lya absorption lines as blends of 3 clouds each with thermal 
width (2kBT/m}iy/'^ ~ 20kms~"^, and a velocity dispersion of llkms""^ (Hu et al. 1995). We 
will also assume that intergalactic material acts only as a "fllter" of external background photons 
emitted from quasars and/or star-forming galaxies, i.e., there are no "local" discrete sources of 
UV radiation within individual clouds. Although the photoionization equilibrium of a H and He 
gas is standard textbook topic, we are applying it in such an unusually wide range of parameters 
(about 7 decades in iVni and nearly 2 decades in J912) that it is worth stating the basic equations 
and approximations used. 
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3.1. Clouds Optically Thin To LyC Radiation 

If we define Ym = iVm/iVH, inei = iVHei/iVne and FHeii = iVHeii/iVne, the steady state 
ionization equilibrium equations can be written as 

1 — ^Hi = Ymlm 

^Hell = ^Hel/Hel , (8) 
1 — ^Hel — ^Hell = iHell/Hell 

where 

Im ^ -^^y (9) 

Here, 

rHi= r dvATT-^amiv) (10) 

is the photoionization rate per hydrogen atom, ckhi the (Case A) recombination rate coefficient to 
all levels of hydrogen, ng the electron number density, and J the background radiation intensity. 
Analogous expressions are valid for /nei and /nell- CoUisional ionizations have been neglected: we 
will discuss in § 3.4 the validity of this approximation. The radiative recombination coefficients 
for hydrogenic ions are taken from Arnaud & Rothenflug (1985), while those for neutral helium 
are taken from Black (1981). The H I and He II photoionization cross-sections are adopted from 
Osterbrock (1989). For the photoionization cross-section of neutral helium we used the following 
fit (to within a few percent up to 1 keV) of the results of Reilman & Manson (1979): 

0.694 X 10"^^ 2 / ^ 

^Hel - ^1.82^^3.23 ' W 

for E > 0.246, where E is the photon energy in units of 100 eV. 

Notice that, while He III recombines for J228,-24 ~ ?^ll,-4, hydrogen is mostly ionized for 
/9i2,-24 ^ 0.1 nn, -4. Thus, in the case of a QSO-dominated background, the typical conditions 
of the Lya clouds assure that most of the gas is completely ionized. In this case we can write 



where 

^ iVHell 
^thin — ,j- 

^Vhi 



r Nmic^m + ??aznC^Heii), (12) 

nHe (1 -MHl)i"HeI ^, ^ f ^HI ^ ^-^r^^ 



thin nn 1 + /Hel(l + /Hell) 12V/HeII. 

assuming a helium to hydrogen cosmic abundance ratio equal to 1/12. 

The solution of the ionization equilibrium given above neglects any kind of radiative transfer 
within the individual clouds, and has been extensively used in previous investigations (Miralda- 
Escude & Ostriker 1990, 1992; Madau 1991, 1992; Giroux & Shapiro 1995). These authors 
have argued that equations (12) and (13) might describe reasonably well the overall transmission 
properties of the universe even when optically thick clouds are present. The reason is that, for 
AGN-like or softer spectra, the ratio rjthin ~ 1-8J912//228 is much greater than unity: this implies 
that clouds which are about to become optically thick in the He II LyC are still optically thin 
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in the H I LyC. To be more specific, let us denote with iVHi(4Ryd) the H I column at which 
He II self-shielding sets in. Clouds with iVni = iVHi(4Ryd) have optical depth of unity at 4Ryd, 
or, equivalently, an He II column iVHell(4Ryd) = crggjj(4Ryd) = 6.35 X lO^^^cm"^. (Similarly, 
self-shielding at 1 Ryd sets in above iVHi(lRyd) = crgj^(lRyd) = 1.59 X 10"^^cm~^.) While our 
estimate of the singly ionized helium to neutral hydrogen ratio from equation (13) is then quite 
accurate for all values of iVni less than iVHi(4Ryd), it is clearly inadequate at larger H I columns. 
In this regime, the absorption of He II LyC photons within the clouds must be taken into account 
in calculating the expected ionization structure. It turns out, however, that the contribution of 
clouds with iVni > iVHi(4Ryd) to the effective optical depth of the universe at 4 Ryd is more or 
less independent of their actual He II column. This is because the exact value of r for r greater 
than a few is irrelevant in equation (5), as most of the photons are removed anyway. For a column 
density distribution dN / dN^n oc iV^j ' like the one given in equation (7), most of the effective 
opacity at frequency v actually comes from absorbers with t{v) ~ 1. 

We have found that, although the optically thin formula correctly estimates the cosmic 
transmission at (or just above) the He II Lyman edge, it breaks down when A ^ 228 A. The 
problem is that, due to the wavelength dependence of LyC absorption, absorbers which are just 
optically thick for these penetrating photons [r(A ^ 228 A) ~ 1] are very optically thick at 
4 Ryd. In these clouds, the ratio iVnell/^Hl is generally much higher than the value estimated 
from equation (13), since He III almost completely recombines (see below). The trasmittance of 
photons above 4 Ryd actually drops to zero for clouds with iVni > iVHi(4Ryd). The net effect is 
that the Universe will be more opaque above 4 Ryd than previously estimated. In the next section 
we will provide a more accurate, albeit still approximate, radiative transfer solution for the actual 
ionization state of individual optically thick absorbers. 

3.2. Clouds Optically Thick to LyC Radiation 

We shall consider a pure H/He semi-infinite slab. The ionization equilibrium equation for 
the neutral hydrogen density at depth R within the slab is 

4:Tr J(iy) _^(^) f 

am[T)nenmi = nm / du am[T^) — ; e ' + 

Jut. hv 

where hv^ = 1 Ryd, 

t{v)= / c/r [nHlCrHl(/^) + nHelC^Hel(i^) + ?^HeIlC^HeIl(i^)], (15) 



dvpaUvf-4^. (14) 

hv 



and the term pabs is the fraction of recombination continuum radiation with emissivity which 
is absorbed within the slab. Analogous expressions can be written for the helium ions. Since the 
recombination emissivity is strongly peaked at threshold, we can approximate the last term in 
equation (14) as 

dvpabs{v)—, ^ Pa6s(l %d)ai25nenHII, (16) 

hv 

where ai2s is the recombination coefficient to the ground level of hydrogen (see Appendix A). For 
simplicity, and because of the steep dependence of the photoionization cross-section with energy, 
we have neglected local helium recombination radiation as a source of hydrogen-ionizing photons. 
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Note, however, that hehum recombination photons wiU be more effective in ionizing hydrogen as 
they escape into the intergalactic space, and are redshifted by the cosmological expansion (cf § 4). 
We shaU take 

1 _ g-r(lRyd) 

P„UiRyd) = i- . (17) 

Equation (17) is strictly vahd only in regions of constant ion and electron density. ^ It is 
possible to show that, for the conditions typically encountered, the product ngnHii is indeed 
fairly constant even throughout absorbers which are very optically thick (log iVni ~ 19) to 
hydrogen LyC radiation. This is not true for helium, however, as He III recombines faster than 
hydrogen, and the radiation intensity at the He II edge is in general smaller than the intensity at 
the H I ionization edge. The He III number density in clouds which are optically thick at 4 Ryd 
is constant up to the He II photosphere, where nHellcilell(4 Ryd)i? 1, then drops sharply by 
more than an order of magnitude as He III starts to recombine. Since most of the reemission and 
absorption of photons above 4 Ryd occurs at the He II photosphere, and r(/ii^ > 4 Ryd) rHell(i^), 
the absorbed fraction of He II recombination radiation can be written as 

=1= 

Pa&.(4Ryd)^i - " , (18) 



where 

■R 



T* = min[ / (irnHeii(?^)c^Heii(4Ryd), 1]. 
Jo 



(19) 



Therefore, while only a fraction r~"^(lRyd) of H I LyC photons will escape from clouds which 
are very optically thick at 1 Ryd, at least 60% of the He II LyC photons produced in individual 
systems will always be able to leak out into the IGM. 

We can now solve for the nni, ?^iiel, and nnell number densities as a function of R. We obtain 

nm _ Ue am .^qx 

~ /;^t//.^aHi(//)e-(-) +neaHl' 

where a}ii{T) = ctm — Pa6s(l Ryd)ai25. Again, similar expressions can be written for the nHel/?^lle 
and njiell/n-Re number ratios. Finally, the H I , He I , and He II column densities can be derived 
by integrating through the slab. 

3.3. Simple Approximations 

An iterative solution to the radiative transfer equation (2) which included a detailed com- 
putation of the ionization structure of individual absorbers would be very time consuming. In 
this section we develop a simple approximation to the exact form of the function rj = He II /H I 
which is capable of reproducing the full absorption properties of intergalactic material. 



^ Assuming a spherical geometry for the clouds would change eqs. (16) and (17), allowing 
slightly more photons to escape for a given optical depth. For r ^ 1, pabs ~ 1 independent of 
geometry. When r ^ 1, the difference between a slab and a sphere can be significant, but in this 
case Pabs ^ 1 anyway, and its exact value is irrelevant in the definition of a. 
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The asimptotic limit for low iVni is readily obtained in the optically thin regime (see eq. 
[13]). The self-shielding of He II LyC radiation becomes important as iVni approaches the value 

iVHi(4Ryd) = 4iVHi(lRyd)ry-i,„. (21) 

At about twice this value He III recombines and the He II column rises sharply, while hydrogen is 
still mostly ionized. The He II to H I column density ratio quickly reaches its maximum value, 

^ nHe 0.37 Thi + neamjT = 1) . . 

nn rieamyT = I) 

where the factor 0.37 (~ 1/e) takes into account the partial self-shielding of neutral hydrogen, Fhi 
is the photoionization rate in the limit of an optically thin cloud, a}ii(T = 1) = ckhi — 0.37ai25ne, 
and we have assumed a negligible fraction of neutral helium. It is straightforward to write down 
a more general expression including the He I component. The contribution of He I to the overall 
opacity of the universe is fully taken into account in our numerical calculations. This is found to 
be negligible except at the very highest redshift, 2; > 5, where the quasar emissivity drops and 
helium starts to recombine (see § 5 below). Finally, at even higher H I columns, iVni > several 
xlO^^ cm~^, both H 11 and He II recombine. The He II to H I ratio becomes a decreasing function 
of H I and finally reaches the asymptotic limit rjthick = of a completely neutral gas. 

We have found that the exact form of //(iVni) can be well approximated by a step-like function, 
as shown in Figure 1 for different values of the background intensity, an assumed temperature of 
T = 2.5 X 10"^ K, and a hydrogen gas density of nn = 2 X 10""^ cm~^ (a larger nn would decrease 
the value of rjmax^ thereby decreasing the cosmological effective opacity above 54.4 eV). The drop 
at large hydrogen columns is well modeled by assuming r] oc l/iVni for iVni ^ 6 X 10"^^ cm~^. A 
similar approximation has been used for the neutral helium component. Note that, in order to 
take into account the fact that both faces of the slab are actually illuminated by the metagalactic 
flux, we have associated the value of at a given i? to a cloud with column density equal to 
2iVHi(-R). Figure 2 shows, again as a function of H I column, the absorption probability through 
an individual cloud, 1 — exp(— r), for various ionizing photon energies. We compare the exact 
value of this probability function with those obtained from the "step-like approximation" and 
the optically thin, rj = const = rjthim limit used in previous investigations. When integrated 
over the actual column density distribution, the step-like approximation gives typical errors of 
a few percent, with a maximum value ~ 20% for photon energies > 200 eV. However, at such 
short wavelengths, the effective optical depth of the universe is < 1, and the percental error in 
computing the speciflc background intensity J (oc e~'^^ff ) is a factor ~ Tg// smaller. 

Since the integration over iVni in equation (5) can be done analytically over ranges of con- 
stant r] and for rj oc l/iVni, the step-like approximation for the ionization structure of individual 
absorbers along the line of sight reduces signiflcantly the computational time. This will allow 
a fast computation of the contribution of recombination radiation to the speciflc intensity (see 
Appendix B). Different possibilities for the nature of the ionizing sources and uncertainties in the 
column and redshift distribution of intervening absorbers can then be quickly explored. 

3.4- Importance of Collisional Effects 

In the relevant range of densities for the Lya clouds and the Lyman-limit systems, nn ~ 
10~'^-10~^ cm~^, collisional ionizations are always negligible for the helium component because 
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of the large ionization potentials. For hydrogen, coUisional ionization can be important when the 
photoionization rate is low, i.e. for low background intensities or in the core of optically thick 
clouds. A related issue concerns the importance of coUisionally excited line emission, notably 
H I Lya , relative to recombination radiation. Since coUisional processes are strong functions 
of the gas temperature, it is important to check that they are indeed negligible in the relevant 
regime. We have therefore solved the equation of thermal equilibrium for a gas slab immersed in 
a background radiation field. In Figure 3 (upper panel) we show, for three different background 
intensities, the average temperature of the cloud as a function of H I column density. For illus- 
trative purposes, and to show the dependence on the assumed gas density, the total hydrogen 
density is taken to be 2 X 10""^ cm~^ for the Lya forest clouds (iVni < 1.59 X 10"^^cm~^), and 
4 X 10~^ cm~^ for the thicker Lyman-limit systems. Heating is provided by photoionization of hy- 
drogen and helium, while recombinations and coUisional excitation of neutral hydrogen dominate 
the cooling term. 

Figure 3 (lower panel) also shows the fractional population of the n = 2 level of atomic 
hydrogen due to recombinations of photoionized H II ions, coUisional excitations from the ground 
state, and recombinations of coUisionally ionized H II ions, for clouds of different H I columns. The 
last contribution is never important, assuring that coUisional ionizations can be neglected in the 
computation of the ionization structure of absorbing clouds. This is also seen in Figure 1, where 
the dotted lines depict the derived He II to H I ratio including coUisional effects. Deviations from 
the values computed in § 3.3 by assuming a purely photoionized isothermal slab are indeed small. 
It is fair to note, however, that a fraction ~ 20 — 30% of H I Lya and two photon decay emissions 
could still be due to radiative de-excitations of coUisionally excited atoms. This contribution has 
not been included in our numerical calculations of the metagalactic flux. 

4. RECOMBINATION RADIATION FROM QUASAR ABSORPTION SYSTEMS 

The attenuation due to the accumulated H I and He II LyC absorption by Lya forest clouds 
and Lyman-limit systems causes a large reduction of the background ionizing flux. However, a 
signiflcant fraction of the absorbed photons is reradiated by the photoionized gas in the ultraviolet 
band. In this section we will show that diffuse radiation from intergalactic material will in fact 
dominate the cosmic ionizing emissivity at characteristic frequencies. Candidates for the source 
of the metagalactic flux includes AGNs (see § 5) and star-forming galaxies, but the precise nature 
of the background is not important in this context. 

The dominant emission from a highly photoionized, dust-free gas which contains, by number, 
92% of hydrogen, 8% of helium, and traces of metals, at a temperature of 2-4x10"^ K, is due to 
radiative recombinations and coUisional excitations of the H I and He II Lya lines at 1216 and 
304 A, respectively. As only Lyman series lines of He II ions can signiflcantly affect the ionization 
state of the universe, and since in the considered range of densities coUisionally excited He II line 
emission is always negligible, we have taken into account only radiation arising from radiative 
recombination processes. Under the typical physical conditions encountered, the hydrogen and 
helium ions in the absorbing medium are in their ground state as recombinations to excited levels 
are followed by radiative decays to the ground level. Among direct recombinations and radiative 
decays, only few processes are important as sources of far-ultraviolet photons. In the following 
we analyze them in turn. 

a) He II processes Direct captures of free electrons to the n^L levels of He II are an ob- 
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vious source of ionizing photons. As the resuhing continuum radiation is strongly peaked at 
frequency Vth = Z'^ /n'^VL, direct captures to the ground level will produce a significant amount of 
He II ionizing photons, while captures to the 2^L level will radiate a Balmer continuum just above 
PL- Within clouds which are optically thick to resonant scattering in the lines of the Lyman series, 
all helium recombinations eventually populate the 2^5" and 2^P levels and are converted into two- 
photon continuum and Lya 304 A photons. In the absence of metals and dust the Lya photons 
will diffuse into the wings and finally escape. Most He II recombination radiation will be emitted 
from clouds with iVni ~ iVHi(4Ryd) ~ 5 x lO^^^cm"^, as these are the dominant absorbers of 
photons above 4 Ryd, and will escape into the intergalactic space without appreciable destruction 
by local H I absorption. 

b) H 1 processes Although only free-bound captures to the ground state of hydrogen are able 
to generate photons above 1 Ryd, we have also computed the contribution of H I Lya and two- 
photon continuum to the near-ultraviolet background. As in the case of He II processes, we have 
taken into account only radiation arising from radiative recombinations. As discussed in § 3.4, 
these account for at least 70-80% of the total hydrogen Lya and two-photon continuum emission. 

c ) Emissivity from Radiative Recombinations To summarize, our calculations will include the 
emission from Lya forest clouds and Lyman-limit systems due to H I and He II Lya recombination 
radiation, two-photon continuum, and free-bound captures to the ground and (for He II only) to 
the 2^ L levels. We have ignored He I recombination radiation because of the smallness of the 
^Hel/^Hell and iVnel/^Hl ratios encountered in typical situations. Note that this will be a 
poor approximation at the very highest redshift, 2; > 5, where the abundance of neutral helium 
increases as the space density of observed quasars drops. 

In general, the specific recombination emissivity associated with a given atomic transition 
can be expressed as 

dz 

er{v,z) = hvf{v)B{z,v)—, (23) 

where 

^ r a^ff r . . ,47rJ(//',z) 

Jo dNmdz a Jz2^^ hv' 

This formula simply states that the rate of recombination photons which escape into the inter- 
galactic space is proportional to the number of photons absorbed per second at any given redshift 
by the relevant atom or ion, times the fraction of recombinations which lead to the radiative 
transition under consideration (the ratio a'^-^-^/a, where a'^^^ is the effective recombination coef- 
ficient, Osterbrock 1989), times the fraction pem of the emitted photons which are able to escape 
from an individual cloud. The fraction of incident photons which are absorbed by a given ion X, 
having optical depth rx, can be written as 

w^Uv') = I drxi^yr^^'''^ ^ [1 - e-^(^')] ^4^. (25) 

The last equality is correct as long as the He II to H I and He I to H I ratios are constant through 
the slab, and it is easily generalized if these ratios vary in a step-like manner instead. In analogy 
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with equation (17) and (18), the fraction of He II LyC photons which can escape from a cloud is 

f ^^^^Sf^ [iVHi<iVHi(4Ryd)]; 
l^^rfeTH^S^ [iVHi>iVHi(4Ryd)]. 

Note that for energies > 54.4 eV we have that r(i^) ~ THell(i^)- Only hydrogen and, above 24.6 eV, 
neutral helium will absorb He II Balmer continuum, two-photon continuum, and Lya photons, 
as well as H I LyC. The function f{v)dv in equation (23) represents the normalized probability 
per recombination that a photon is emitted in the interval {v, v + d.v). We stress that the 
recombination emissivity from intervening material is largely independent on the (poorly known) 
sizes and gas densities of the absorbers, and depends only weakly on the temperature of the 
emitting media. 

d.) Recombination Continuum The rate of direct recombinations to the n^L level can be com- 
puted using the Milne relation (Osterbrock 1989). This, together with equation (23) for hydrogen- 
like ions yields 

mo^n^L = % ilT^)''' ^.VHi(./.a)e-'^(^--)/'^^ (27) 

where all symbols have their usual meaning or have been defined in the text. Direct captures 
to the 2^L level of He II result in a radiation continuum with specific volume emissivity equal 
to 2(iVHelll/^llll) ~ 2(iVHe/^ll) ~ 1 / 6 of the emissivity arising from free-bound captures to the 
ground state of H I . 

t) Lya Line Radiation Inserting equation (23) into equation (2), and taking f{ry) = 8{v — Vc(), 
where Va is the Lya frequency and 8{x) is the usual delta function, the specific intensity associated 
with Lya line radiation can be written as {v < Va) 

Ja{vo,Zo) = — / dzv8{v-v^ y._/^ ^(z,//^)e-"-//('-°'^°'^) 
47r I ^ I i H" ^ ) 

'° 2 (28) 
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where Zem denotes the redshift of the Lya line emitters, 1 -|- Zem = (1 + Zo)valvo- The effective 
recombination coefficient to the 2^P level has been assumed to be independent of column density, 
since resonant absorption of Ly/3 photons populates the 3^P level, which subsequently decays via 
Ha and two-photon emission. The population of the 2^P level arising from downward cascading 
of n > 3 levels (following the resonant absorption of higher order Lyman lines) can be neglected 
(see Pengelly 1964 and Appendix A). 

f) Two-Photon Continuum The 2^5" level is populated by direct recombinations and by down- 
ward cascades following recombinations to higher levels. The function gBiiv) = hv f{v) for the 
two-photon continuum decay of neutral hydrogen is tabulated in Osterbrock (1989). The emis- 
sivity due to the two-photon continuum decay of the 2^5" level of He II can be computed using 
the relation 

g^eiiiv) = gmivl^). (29) 
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The effective recombination coefficient to the 2^5" level strongly depends on the Ly/3 line center 
opacity, r(Ly/3). As this increases from to oo, ^225 increases by a factor 2.6. However, for 
clouds with r > 1, which are the major contributors to the recombination emissivity, r(Ly/3) > 
10^. We have therefore assumed that the population of the 2^5" is well described by Case B 
recombination for all values of iVni (see Appendix A for details). 

5. NUMERICAL RESULTS 

The formalism we have developed so far, including the effects of the reprocessing of UV 
photons by a clumpy intergalactic medium, can be applied to compute the metagalactic flux 
arising from any discrete source of photoionizing radiation, the most plausible candidates being 
AGNs and young, star-forming galaxies. The last contribution is the most difficult to assess and 
also the most intriguing, as a determination of the UV background intensity and/or spectrum 
at high redshift might constrain the history of cosmic metal production and galaxy formation 
(Bechtold et al. 1987; Miralda-Escude & Ostriker 1990; Songaila, Cowie, & Lilly 1990; Madau & 
ShuU 1995). The intensity and spectrum of the ultraviolet background arising from massive stars 
in metal-producing galaxies will be the subject of a companion paper (Madau & Haardt 1995). 
Below, we will focus our attention on quasar sources. 

a) Quasar Emissivity as a Function of Redshift According to recent surveys (see, e.g., the one 
by Boyle, Shanks, & Peterson 1988, with about 400 ultraviolet-excess quasars), a pure luminosity 
evolution model, in which QSOs statistically conserve their number since 2; ~ 2, describes well 
the global properties of the quasar population at low redshift. Beyond 2; ~ 2, the comoving space 
density of quasars appears to stay constant up to a redshift of ~ 3 (see, e.g., Koo & Kron 1988; 
Boyle, Jones, & Shanks 1991). The decline in the quasar counts at 2; > 3 has been the subject 
of a long-standing debate, and it is still controversial. It has been argued by Boyle (1991) and 
Irwin et al. (1991) that the space density of bright {Mb < —26) quasars at 2; ~ 4.3 is similar 
to that at 2; ~ 2. On the other hand. Warren, Hewett, & Osmer (1994) find strong evidence 
for positive evolution of the quasar luminosity function (LF) beyond 2; ~ 2. According to these 
authors, the evolution ceases at z ~ 3.3, with a marked decline by a factor > 3 in space density 
at larger redshift, as first suggested by Osmer (1982). A decline in the LF is also seen at the 
low luminosities {Mb ^ —26) probed by the Schmidt et al. (1991) grism survey. We remark that 
these estimates are quite uncertain because of corrections for incompleteness and the adopted 
continuum emission spectrum, and because of the possible bias introduced by dust obscuration 
arising from intervening systems (Ostriker & Heisler 1984; Wright 1990; Fall & Pei 1993). 

We shall estimate the integrated volume emissivity from QSOs following Pei (1995), who has 
recently derived an analytical model which fits well the empirical LF estimated by Hartwick & 
Shade (1990) and Warren et al. (1994). ^ The best-fit to the QSO LF is the standard double 
power-law: 

<f>{L,z) 
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^ Hartwick & Shade (1990) have compiled 15 major optical surveys with more than 1000 
quasars covering moderate redshifts 0.1 ^ 2; < 3.3. Most recently. Warren et al. (1994) have 
combined their own and four other surveys to form a sample with contains about 200 quasars 
covering the range of high redshift 2 < 2; < 4.5. 
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The entire LF shifts along the luminosity axis as the position of the break evolves with redshift: 

L^{z) = L^{0){1 + z)"-^ exp[-z{z - 2z^)/2al], (31) 

where a power-law spectral distribution for the typical quasar UV spectrum has been assumed, 
f{ry) oc and the dependence on the spectral index a is explicitly shown. According to Pei 

(1995), this Gaussian form can fit reasonably well the observational data in the entire range 
< 2; < 4.5; the evolution of QSOs reaches a maximum at 2; ~ 2.8 and declines at higher redshits. 
In a cosmology with h^o = 1 and qo = 0.1, the LF fitting parameters values are /3i = 1.83, 
/32 = 3.7, 2;* = 2.77, cr* = 0.91, and log (^*/Gpc~^) = 2.37; the B-magnitude at the break is 
Mb,* = — 23. At the present epoch this LF is in good agreement with estimates of the Seyfert 
LF by Cheng et al. (1985). The proper volume emissivity can be written as 

eQ(iy, z) = eqivB, 0)(1 + z)«+2(/////b)"" exp[-z(z - 2z,)l2all (32) 

where eq^VB^^) is the extrapolated z = Q emissivity at the reference frequency vb = c/4400A: 

(p{L, 0)LdL ~ 6.4 X 10^^ ergs Gpc"^ s"^ Hz"^ . (33) 

—'min 

The LF of Seyfert galaxies matches remarkably well that of optically selected QSOs at Mb = —23, 
and does not show clear evidence of leveling off down to Mb ~ —18.5 (Cheng et al. 1985), 
corresponding to Lmin ~ 1.5 X 10~^-L*(0). For the given shape of the LF, a value of Lmin 4 times 
smaller would increase the emissivity by only 20%. 

b) Quasar Spectral Energy Distribution One of the biggest uncertainties in these calculations 
has always been the detailed shape of the UV continuum spectrum of quasars. Recently, however, 
the results of several spectroscopic surveys have significantly advanced our knowledge of the QSO 
emission properties. In particular, a very high signal-to-noise composite spectrum of the rest- 
frame ultraviolet and optical region of about 700 high luminosity, intermediate redshift quasars 
(part of the Large Bright Quasar Survey) has been constructed by Francis et al. (1991). Force- 
fitting a power-law to the constituent spectra yields a median spectral index of ckq ~ 0.3 in the 
range 5000-1500 A. In the UV, longward of the Lya emission line, the survey of Sargent, Steidel, 
& Boksenberg (1989) yields a softer slope, ajjv = 0.78 it 0.27 (see also O'Brien, Gondhalekar, & 
Wilson 1988). Finally, the Tytler et al. (1995b) Hubble Space Telescope FOG snapshot survey (82 
QSOs at 2; ~ 3) has provided, for the first time, statistical information on the intrinsic ionizing 
photon distribution of quasar. After accounting for intervening absorption, Tytler et al. find an 
even softer slope, ai = 1.4 it 0.25, in the wavelength range 1300-330 A. 

We shall therefore adopt the following model for the "typical" quasar spectral energy distri- 
bution: 

( v-^-^ (2500 < A < 4400 A); 
/(//) oc i (1216 < A < 2500 A); (34) 

[ (A<1216A), 

where the different slopes have been continuosly matched. In the far-UV, this model is similar to 
the "medium" QSO spectrum used in Bechtold et al. (1987), and "Model QS2" of Miralda-Escude 
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& Ostriker (1990); its extrapolation to 2 keV gives ctox — 1-4, in agreement with the value for 
QSOs measured by the Einstein Observatory (Zamorani et al. 1981). The soft (0.2-2 keV) X-ray 
spectral index is consistent with ROSAT observations of QSOs (e.g., Laor et al. 1994), but is 
significantly steeper than suggested by earlier Einstein and EXOSAT data (Wilkes & Elvis 1987; 
Comastri et al. 1992). Note that we have ignored luminosity and redshift dependent effects in 
QSO UV/soft X-ray emission spectra (O'Brien et al. 1988; Bechtold et al. 1994), and that the 
adopted X-ray spectral slope is too steep to fit the diffuse X-ray background below 40 keV (see, 
e.g., Madau, Ghisellini, & Fabian 1994 and references therein). Also, no attempt has been made 
to include the soft X-ray excess emission below 0.2 keV which has been detected in many AGNs 
(Turner & Pounds 1989) and is well modeled by the tail of a thermal spectrum with kT ~ 40-80 
eV (Comastri et al. 1992). It has been suggested (see Sciama 1994) that the presence of this 
thermal, far-UV bump might help to account for the large abundances of N V and O VI observed 
in some Lyman-limit systems. 

In Figure 4 we plot the (proper) volume emissivity at the H I , He I , and He II Lyman edges 
due to QSOs as a function of redshift. Note that the "A'-correction" term in equation (32) has 
been matched to take into account the varying spectral index a of equation (34). In the redshift 
range < 2; < 1, e(912A,2;) matches to within 30% the volume emissivity adopted in Madau 
(1992) from the Boyle (1991) quasar LF. At z ^ 1.5 (3) the former is about 2 (2.5) larger than 
the latter, then falls like a gaussian beyond 2; ~ 4 in contrast to the constant comoving emissivity 
scenario (Boyle 1991; Madau 1992). 

c ) Spectrum of the Ionizing Background The formal solution of the radiative transfer equation 
(2) is not an explicit form of the specific background intensity, as its right-handed term implicitly 
contains J in the emissivity term, e = eg -\-€r{J)^ and in the effective optical depth Teff{J). Phys- 
ically, this simply means that the metagalactic flux depends on the ionization state of intervening 
clouds which is in turn determined by the background fleld itself. In order to solve for J, an iter- 
ative scheme must be implemented (see Appendix C). Figure 5 shows the background radiation 
fleld produced by the observed QSOs in a go =0.1 cosmology, as a function of redshift from the 
present epoch to z = 5.5. The spectrum from 5 to 5000 A results from the numerical integration 
of equation (2), and is self-consistently flltered through a clumpy IGM. For comparison we have 
also plotted the same spectrum in the limit of a perfectly transparent (te// = 0) universe, and 
a purely absorbing intervening medium. As already discussed by a number of authors (Bechtold 
et al. 1987; Miralda-Escude & Ostriker 1990, 1992; Madau 1991, 1992; Meiksin & Madau 1993; 
Zuo & Phinney 1993), the absorption by hydrogen and helium in intervening material produces 
a very strong effect in the spectrum of the ultraviolet metagalactic flux. This will be generally 
characterized by two breaks at the H I and He II ionization edges, the amplitudes of which gener- 
ally increase with redshift. The background intensity above 4Ryd is, in particular, a non-linear 
function of the number of photons emitted by QSOs at these frequencies: a relative small de- 
crease in this number causes a large increase in the size of the He II break. The importance of the 
reprocessing of ionizing photons by intergalactic material is clearly seen in Figure 5. While red- 
shifted He II Lya and two-photon continuum emission dominate the metagalactic flux from ~ 1 
to 3 Ryd, direct recombinations to the ground level of hydrogen and singly ionized helium provide 
a substantial number of LyC ionizing photons close to the ionization thresholds. Because of the 
larger cosmic opacity at early epochs, the amplitude of the jump in the background intensity at 
1216 A and 304 A due to H l and He II Lya emission increases with redshift from 2; = to 2; 4, 



to decrease again at 2; > 4 as the QSO emissivity drops, and the background flux becomes dom- 
inated by more local, less attenuated sources. At low redshift, Lya line recombination radiation 
forms an extended hump in the background spectrum due to the integrated contribution of QSO 
absorbers along the line of sight. The diffuse Lya flux becomes more peaked at early epochs, 
again because of the rapid decrease in the quasar ionizing emissivity above 2; ~ 4. Also note how, 
at 2; > 5, helium starts to recombine: this produces a small absorption feature in the spectrum 
of the metagalactic flux at 504 A. 

d) Photoionization Rates The net effect of recombination emission from QSO absorbers on 
the ionization state of the universe is best depicted in Figure 6, where the total photoionization 
rates rHi,rHel, and FHell are plotted as a function of redshift and compared to the case of a 
purely absorbing IGM. ^ Relative to the latter, our model calculations including riemission 
produce a value of Fhi about 1.5 times higher, nearly independent of redshift. Typically, less 
than 1/3 of the IGM diffuse contribution comes from optically thick Lyman-limit systems. The 
effect of recombination radiation makes Fnell 1-6 times larger at 2; = 1, and 2.2 times larger at 
z = 3. How do these rates compare with observational constraints? Diffuse Ha surface brightness 
observations of the Haynes-Giovanelli intergalactic cloud by Vogel et al. (1995) have recently 
yielded an upper limit to the local metagalactic photoionization rate of 3 X 10""^^ . The total 
QSO+IGM contribution at 2; = from Figure 6 is about eight times lower, Fhi — 4 X 10""^ . 
Because of the larger QSO emissivity adopted here, and the inclusion of recombination radiation 
from intergalactic material, this value is about two times higher than the one obtained by Madau 
(1992). The amount of diffuse ionizing radiation predicted by our model at the present epoch 
is consistent with the value derived by Kulkarni & Fall (1993) by applying the proximity-effect 
technique towards QSOs at {z) 0.5. Observations at 21 cm of the fairly abrupt truncation in 
the H I distribution at the edge of nearby spiral galaxies have also been used to infer values of the 
metagalactic photoionzing flux. This technique is somewhat model-dependent, and yields rates 
in the range 10"^^ - 3 X 10"" s'^ (Maloney 1993; CorbeUi k Salpeter 1993; Dove k ShuU 1994), 
again consistent with our model calculations. 

At high redshift, the proximity effect, which is the measured decrease in the number of Lya - 
absorbing clouds induced by the UV radiation fleld of a QSO in its neighborhood, provides an 
estimate of the intensity of the metagalactic flux at the hydrogen Lyman edge. Determinations 
of J912 in recent years include the following: log J912 ~ —20.5 at 1.6 < z < 4.1 (Bechtold 1994), 
between -22 and -21.5 at z ^ 4.2 (Williger et al. 1994), log J912 = -21 ± 0.5 at 1.7 < z < 3.8 
(Bajtlik, Duncan, k Ostriker 1988; Lu, Wolfe, k Turnshek 1991), and log J912 = -21.3 ± 0.2 at 
1.7 < z < 3.8 (for a velocity offset of 1400 km s~"^ between the emission lines in the broad line 
region and the QSO rest-frame, see Espey 1993). If J,^ oc , the corresponding photoionization 
rate equals 3 X 10""^^ J9i2,-2i . Our QSO-dominated background model generates Fhi ~ 
1.6 X 10""^^ at z = 2.5, and a value three times lower at z = A. We conclude that, within the 
uncertainties, the observed QSOs appear able to provide the number of ionizing photons required 
by the proximity effect at z ^ A. 

^ A Gaussian flt, F = A{1 -\- 2;)^exp[— (2; — Zc)"^ j S]^ yields for the H I photoionization rate 
plotted in Figure 6 A = 6.7 X 10"" sec'^ B = 0.73, z^ = 2.30, and S = 1.90. This flt can 
be applied over the range < 2; < 5 with a maximum error lower than 10%. A similar fltting 
formula to 10% accuracy in the same redshift range gives for the He II rate A = 6.7 X 10""^^ sec~"^, 
B = 0.72, Zc = 2.35, and S = 1.98. 
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e) Integrated H I Lya Emission at z = from Quasar Absorption Systems We now address 
the long-standing question of whether redshifted 1216 A H I Lya emission from intergalactic gas 
might be detected in the diffuse far-UV background in the range 1300-2000 A. The intensity 
of the ionizing diffuse flux from quasar sources imphes that the redshift-smeared background 
at 2; = due to H I Lya recombination radiation from photoionized Lya clouds and Lyman- 
limit systems, about 3photcm~^ st~^ at 1500 A, is at a level of less than flve percent of 
current observational limits (see Bowyer 1991 for a review). Therefore, in a universe photoionized 
by QSOs, redshifted Lya recombination emission from intergalactic material can be ruled out as 
a signiflcant source of far-UV radiation. 

f) Intensity Fluctuations of the Ionizing Background The attenuation by QSO absorption- 
line systems reduces signiflcantly the total number of primary sources involved in the production 
of the metagalactic flux. The qualitative behaviour can be readily obtained by expanding the 
effective optical depth in redshift around Zg. From equation (B9) we obtain 

Are//(912 A, z) ~ [0.12(1 + zf-^"^ + 0.25(1 + zf-^^] Az (35) 

at the hydrogen Lyman edge, where the two terms on the right-hand side represent the con- 
tribution of the Lya clouds and Lyman-limit systems to the effective opacity. At 2; = 3, 
Are//(912A) ~ 1 for Az ~ 0.17; in a low qo universe, this corresponds to a (proper) absorp- 
tion length of A£ ~ Mpc. The ionizing radiation fleld is thus largely local at 1 Ryd, as 
in the "attenuation volume" 47r(A^)^/3 the QSO LF (eq. [30]) predicts only ~ 4 objects with 
L > L^z = 3). 

As argued by Zuo (1992) and Fardal & ShuU (1993), this fact might have an interesting 
consequence, namely intergalactic absorption should greatly enhance the intensity and spectrum 
fluctuations of the ultraviolet background at high redshift. The importance of this effect re- 
sides in the possibility that the study of the influence of a fluctuating metagalactic flux upon 
Lya absorption lines (e.g., the production of clearings and voids in the forest along the line of 
sight) could shed some light about the nature of the main contributors of ionizing background 
photons (Kovner & Rees 1989). Here we want to stress that, as recombination radiation from 
the numerous (relative to QSOs) Lya clouds and Lyman-limit systems contributes ~ 30% of the 
total photoionization rate at 2; ~ 3, the UV background will be much smoother in space than 
previously considered. In particular, the two-point correlation function induced in the Lya forest 
by intensity fluctuations will be reduced by a factor of about 2 (Fardal & ShuU 1993). 

H I Gunn-Peterson Test In the presence of uniform material of density n(z) along the path 
to a distant object, any resonance line with wavelength A will produce a Gunn-Peterson (1965) 
absorption trough with optical depth Tap{z){= rre'^ /m(,c)II(z)~^ Xfn(z), where z is the redshift 
at which the observed radiation passes through the resonance, II(z) is the Hubble constant, 
/ is the oscillator strength of the transition, and all other symbols have their usual meaning. 
The H I optical depth is related to the photoionization rate through the condition of ionization 
equilibrium 

T^l = 0.18h-,%-'-''{i}ohlof{l + zf{l + 2qozr'/'T^l_,,, (36) 

where T4 is the temperature of the ionized diffuse intercloud medium, its mass density param- 
eter, and rHi,-i2 is the ionization rate in units of 10""^^ . The strongest limit on the amount 
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of diffuse intergalactic neutral hydrogen at 2; ~ 3 is provided by a direct estimate of the quasar 
Q2126-258 flux in regions of the spectrum where hues are absent (GiaUongo, Cristiani, & Trevese 

1992). The derived la upper bound is r^Sp < 0.04. Inserting this hmit into equation (36) with 

3/2 

T4 = 2, z = 3, qo = 0.1, and rHi,-i2 = 1-6 yields Qjjh^l^ < 0.027. Nucleosynthesis constraints on 
the cosmological baryon density require O^^so > 0.04 (Walker et al. 1991), and would be consis- 
tent with the limit on Q]j above if less than 70% of the baryons were in un uniformly distributed 
medium, the rest having collapsed into galaxies and discrete absorption systems. Two possible 
problems must be noted here: 

1 . If the diffuse IGM is photoionized by quasars sources only, the Gunn-Peterson test becomes 
more and more difficult to satisfy at earlier epochs, both because of the strong dependence of tgp 
with redshift and of the rapid drop in the QSO emissivity beyond 2; ~ 4. As depicted in Figure 
6, the background ionization rate due to QSOs falls by more than a factor of 5 from z ^ 3 to 
z 4.3, while the cosmological term in equation (36) increases by 4. Therefore we would expect 
the Gunn-Peterson optical depth to increase over 20 times. However, an analysis of the spectrum 
of the QSO 1202-0725 by GiaUongo et al. (1994) yields tgp < 0.05 at z ~ 4.3. 

2. In the gravitational instability scenario for the formation of cosmic structures, the as- 
sumption of a uniform intercloud medium has no theoretical justiflcation. As recently argued 
by Reisenegger & Miralda-Escude (1995), a photoionized IGM must be inhomogeneous on scales 
larger than the Jeans length, and the Gunn-Peterson opacity should fluctuate. Moreover, as the 
flUing factor of underdense regions or voids is large, the median tqp should be signiflcantly lower 
than the value expected for a uniform medium with Q]j = Qb (Reisenegger & Miralda-Escude 
1995). 

h) He II Gunn-Peterson Test The recent iT^T observations of quasars Q0302-003 at z = 3.29 
(Jakobsen et al. 1994) and PKS1935-692 at z = 3.19 (Tytler et al. 1995b), and the ASTRO-2 
HUT spectrum of quasar HS1700-|-6416 at z = 2.72 (Davidsen et al. 1995), show the presence 
of absorption troughs at the redshifted wavelength of He II 304 A. It is still unclear whether 
line blanketing from discrete Lya forest clouds alone can produce the observed flux reduction, 
or whether a signiflcant contribution from absorption in a smoothly distributed IGM is needed 
instead (see, e.g., Songaila, Hu, & Cowie 1995). In the latter case, the He II Gunn-Peterson 
test could potentially provide stringent constraints on the amplitude and spectral shape of the 
UV ionizing background (Meiksin & Madau 1993; Miralda-Escude 1993; Madau & Meiksin 1994; 
Giroux et al. 1995). The He II and H I optical depths can be related as 

^GP _ 'Ithtn /„„x 

^Hi - 4 • y-^') 

^GP ^ 

In Figure 7 we plot as a function of redshift the He II /H I ratio for an optically thin gas in 
ionization equilibrium with our model ultraviolet background. The predicted value increases 
from 25 at 2; = to ~ 45 at 2; = 2.5, to decrease again below 30 at 2; > 4.5. This is because 
at low redshifts the universe is nearly transparent, while at 2; > 4.5 the decline in the quasar 
counts makes the metagalactic flux dominated by more local, unattenuated sources. Note how 
recombination radiation from intergalactic material tends to decrease the expected value of rjthin 
for z > 1. 

A background dominated by quasar sources could then produce a He II Gunn-Peterson optical 
depth Tqp^ ~ 0.5 at 2; ~ 3 and still be consistent with the lack (within the uncertainties) of any 
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detectable T^Sp. We conclude that the He II 304 A absorption decrement observed by Jakobsen 
et al. (1994), Tytler et al. (1995b), and Davidsen et al. (1995) should be partly due to uniformly 
distributed He II , unless most of the baryons have been removed from the IGM, possibly into the 
Lya forest. 

i) He II Lya Clouds With the launch of HUT on ASTRO-2 it has become possible to detect 
and identify He II Lya lines from Lya forest clouds (Davidsen et al. 1995). These could provide, 
together with the corresponding H I Lya lines, constraints on the spectrum of the metagalactic 
flux. A curve of growth analysis would enable us to estimate the parameter Tjnii^n 

and thus to 

test the QSO-dominated background scenario. Note that, because of the relatively soft diffuse 
radiation fleld produced by QSOs, He II lines associated with H I lines which are still in the linear 
part of the curve of growth are expected to be heavily saturated, even more so if massive stars in 
star-forming galaxies provide a signiflcant contribution to the background flux. 

i) He II Lya Line Blanketing At wavelengths shortward of the H I and He II Lya lines, the 
quasar's continuum intensity is attenuated by the combined blanketing effect of many absorption 
lines arising from intervening discrete systems. As mentioned above, the relative strength of 
He II versus H I continuum and line opacity depends on the intensity of the metagalactic flux at 
4 Ryd compared to 1 Ryd. The recent detection of He II Lya absorption in the spectra of three 
high- 2; quasars can thus be used as a diagnostic tool of the ambient physical conditions at early 
epochs. The effective H I line-blanketing optical depth due to Poisson-distributed absorbers is 
given by 

1 _|_ 7 fWmax pp. AT 

(Paresce, McKee, & Bowyer 1980), where N/ dWdz is the observed rest equivalent width dis- 
tribution. The He II effective optical depth may be obtained from equation (38) using a curve of 
growth analysis and a distribution of Doppler parameters b for H I and He II . From the numerical 
study of Madau & Meiksin (1994), we derive for thermally broadened clouds {b = const ~ 30 km) 
with W > W^^n = 0.0025 A (iVm > 5 x 10^^ cm'^) 

rfffiz = 3.3) ^ 0.35ry°,^^, (39) 

a flt to 3% accuracy over the range 10 < rjthin < 100. As our model background generates 
Vthin ~ 45 (corresponding to a "softness" parameter Sl = J912I J22» ~ O-QVthin ~ 30) at z ^ 3.3, 
we conclude that a He II Lya line blanketing opacity greater than unity is to be expected in 
a universe photoionized by QSOs (Miralda-Escude 1993; Madau & Meiksin 1994; Giroux et al. 
1995). Note that, as the blanketing opacity is dominated by those lines which lie at the transition 
between the linear and the flat part of the curve of growth, T^jj^ is in fact sensitive to the 
undetermined H I distribution at very low H I equivalent widths, and increases if the hydrogen 
and helium absorption lines have the same width (turbulent broadening) instead (Madau & 
Meiksin 1994; Cowie et al. 1995). 

k) The Redshift Evolution of the Lya Clouds As mentioned in § 2.2, the number of Lya forest 
lines seen by the HST at 2; ~ is larger than predicted by extrapolating from the high- 2; regime 
(Bahcall et al. 1991; Morris et al. 1991). It has been suggested by Ikeuchi & Turner (1991) that 
this excess might be related to the sharp drop in the ionizing photon flux expected for z < 2 in 
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a QSO-dominated UV background. We point out that a similar effect should be observable at 
2; > 4. As shown in Figure 6, the photoionization rate Fhi from QSOs is estimated to decrease 
by a factor of 6 between z = 2.5 and z = 4.3. If there is no compensating change in the evolution 
of cloud properties, this drop should produce an increase in the neutral hydrogen column and 
therefore in the number of Lya clouds observed per unit redshift. A strong increase in the rate 
of incidence of Lya clouds at large redshift has been claimed by Williger et al. (1994), who find 
dN/dz oc (1 + zf-^ in the range 3.7 < z < 4.2. 

I) Photoionization Models for Metal-Line Absorption Systems It has been long recognized that 
it is possible to model the ionization structure of the heavy-element QSO absorption systems and 
discriminate between different background spectral shapes by looking at the line strengths of 
various elements (Chaffee et al. 1986; Bergeron & Stasinska 1986; Steidel & Sargent 1989). Mea- 
surements of the column densities of H I and several ionized species of C, N, and O in three 
optically thin systems by Vogel & Reimers (1993) have been used in conjuction with photoion- 
ization models to rule out star-forming galaxies as the main source of the UV background. The 
dominance of high-ionization species can be best explained by an AGN-type spectral energy distri- 
bution with softness Sl < 5. However, when combined with the presence of strong He I resonance 
lines (Reimers & Vogel 1993), these observations appear incompatible with photoionization equi- 
librium models (Giroux, Sutherland, & ShuU 1994). Using Si II , Si IV , and C IV lines, Songaila 
et al. (1995) have recently estimated Sl to be in the range between 20 and 100 at z ^ 3.2. 

Here, we only want to point out the inadequacy of photoionization models which assume 
a background power-law ionizing spectrum. As clearly shown in Figure 5, the reprocessing of 
LyC photons by intergalactic material introduces bumps, wiggles, and breaks at characteristic 
wavelengths in the spectrum of the metagalactic flux, the presence of which must be taken into 
account when modeling metal systems. The integrated background due to He II Lya line emission 
at 40.8 eV is clearly going to affect the abundance of, e.g., O II ions. We flnd that the O II /C III 
ratio is particularly sensitive to the spectrum of the metagalactic flux at 35 eV compared to 
48 eV, and hence to He II Lya recombination radiation (see Fig. 6). The absorption line O II 833 
(together with C III 977) has been recently identifled by Vogel & Reimers (1995) in the spectrum 
of the quasar HS1700-F6416. 

m) A Case with qo = 0.5 For 2; < 1, optical surveys are complete to a sufficiently faint 
magnitude level to detect the QSOs that contribute the bulk of the ionizing background, and 
the computed metagalactic flux does not depend on the assumed cosmology. However, at larger 
z this will no longer be true, as more and more quasars fall below the survey limit. It is then 
necessary to correct for the missing QSOs: the selection function one adopts depends on the value 
of qo. In the attempt to bracket the uncertainties on J{v,z), we have run a go = 0.5 background 
model using the quasar LF fltting parameters given in Pei (1995) for this cosmology. In a flat 
universe, the usual Gaussian flt, F = A{1 -\- z)^ exp[— (2; — Zc)"^ j S] (good to 10% over the relevant 
redshift range) yields for the H I photoionization rate A = 6.7 X 10""^^ sec~"^, B = 0.43, Zc = 2.30, 
and S = 1.95. Similarly, the parameters for the He II rate are A = 6.3 X 10""^^ sec~"^, B = 0.51, 
Zc = 2.3, and S = 2.35. Relative to the go = 0.1 case (see footnote [3]), the flat universe model 
provides about 30% less hydrogen- and helium-ionizing photons above z = 2. 

n) A Case with Low Lntergalactic Absorption As already mentioned in § 2.2, a steepening in the 
iVni power-law distribution of the Lya clouds at high columns, suggested by several observations. 
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may result in a smaller effective cosmic opacity and thus in an increased ionizing metagalactic 
flux (Meiksin & Madau 1993). In order to quantify this eff'ect, we have run a qo = 0.1 background 
model using the same column density distribution of intervening absorbers given in equation (7), 
but assuming no forest clouds between 10"^^ and 1.59 X 10"^^cm~^. The Gaussian flt to the new 
photoionization rates yields A = 7.0 X 10""^^ sec~"^ , B = 0.8, Zc = 2.22, and S = 1.95 for H I , 
and A = 1.4 X 10"^^ sec'^ B = 0.69, = 2.32, and S = 2.00 for He ii . Relative to the 
standard qo = 0.1 case (see footnote [3]), the low attenuation model provides about 30% (20%) 
more hydrogen-ionizing photons at 2; = (2). However, the eff'ect on the He II -ionizing flux is 
much larger, as the lack of Lya forest clouds above 10"^^ cm~^ signiflcantly decreases the cosmic 
opacity above 54.4 eV, thereby increasing Fnell by about a factor of 2. Such a model would be 
characterized by a smaller He II over H I ratio, and may therefore have difliculties in accounting 
for a signiflcant absorption trough at 304 A (Jakobsen et al. 1994). 

6. SUMMARY 

In this paper we have shown that, beside the attenuated direct radiation emitted from quasars 
and/or star-forming galaxies, the ionizing metagalactic flux will also have a signiflcant diffuse 
component due to recombinations within the clumpy intergalactic gas. In a photoionized universe, 
the following processes will contribute to the diffuse hydrogen-ionizing background fleld from 
Lya forest clouds and Lyman-limit systems: 

1. Recombinations to the ground state of H I and He II . 

2. He II Lya emission (2^ P l^S) at 40.8 eV. 

3. He II two-photon (2^5" — > 1^5") continuum emission. 

4. He II Balmer continuum emission at > 13.6 eV. 

In Case A, a fraction ~ 0.7 of the excited-state recombinations will populate the 2^P state, while 
the remainder end up in the 2^5" state. In the presence of a power-law emissivity at a level 
currently envisaged from QSOs, the He II present in QSO absorption-line systems will efliciently 
degrade soft X-ray, He-ionizing photons into ultraviolet, H-ionizing photons. The addition of 
He I recombination radiation would be an improvement of this work, although this should only 
change our results at the very highest, z ^ 5 redshift considered in this paper. We have modeled 
in detail the ionization state of intervening clouds and the propagation of ultraviolet radiation 
through intergalactic space, and shown that the universe is more opaque above 4 Ryd than previ- 
ously estimated. We have also reassessed the contribution of the QSOs observed in optical surveys 
to the UV extragalactic background, and found that the attenuated direct quasar emission plus 
the recombination radiation from the intergalactic gas can provide enough hydrogen-ionizing pho- 
tons to satisfy the proximity effect at large redshift. The reprocessing of LyC background photons 
by the intergalactic medium has been shown to signiflcantly affect the intensity, spectrum, and 
fluctuations properties of the metagalactic flux. 

The recent discovery by Cowie et al. (1995) and Tytler et al. (1995b) of metals in the 
Lya clouds has shown that the IGM at high redshift is contaminated by the product of stars. In 
the context of our results, these measurements raise two important issues, which will be discussed 
in detail in a companion paper (Madau & Haardt 1995): 

(1) As Lya clouds have large flUing factors and contain a signiflcant fraction of the baryons 
in the universe, they might trace the bulk of star formation occurring in galaxies at high redshift. 



- zo - 



Madau & ShuU (1995) have shown that the background intensity of photons at 1 Ryd that ac- 
companies the production of metals in the Lya clouds may be significant, comparable with the 
estimated contribution from QSOs (unless a large fraction of the UV radiation emitted from stars 
cannot escape into the IGM, or if most of the metals observed at z ^ 3 were produced at much 
earlier epochs). If this is true, the signature of a stellar spectrum should appear in several places, 
notably the He II Gunn-Peterson effect and the distribution of metal ionization states of the Lya 
forest and Lyman-limit absorbers. 

(2) The new data could also imply the presence of a widespread distribution of dust in the 
IGM. Even a small amount of dust in the Lya forest and Lyman-limit systems might provide a 
significant source of absorption of ionizing radiation in the universe. Further studies of this effect 
would be an improvement of the present work. 
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useful discussions, and N. Panagia for clarifying various issues concerning the physics of radiative 
transfer. This research was partially supported by the Director's Discretionary Research Fund at 
STScI. FH also acknowledges financial support by the Italian MURST and by the Swedish NFR 
Council. 



APPENDIX A 

RECOMBINATION RATE COEFFICIENTS 
In this Appendix we provide some useful fitting formulae to the hydrogen recombination rate 
coefficients a„2^ (in units of cm^ ^~^) needed for our computations. The recombination coefficient 
to the ground level of hydrogen can be written as 



«i^s^l-58xl0-"x L (Al) 



-13 J T,-'-'' (T4 < 1); 

21-0.51-O.lOln T4 ^rj^^ y -^ -j 

This expression fits the results of Burgess (1964) within 3% for 0.1 ^ T4 < 10. We have also used 

the following fitting formulae (within 2% for 0.1 < T4 < 10) to the Case B results of Pengelly 
(1964): 

{Y^-0.92-0.101nT4 /rp ^ -\ \ 

92 nnfi, T " (A2) 

2^ — 0.92 — 0.06 m J4 (2^>1) 



au' 



d 



I Y^-0.65-0.021n T4 /rp < 1 \ 

al(f ~ 8.37 X 10-1^ x\ n 07, t ~ (A3) 

2^ b I — 0.72 — 0.07m J4 ^rp > ]^) 

APPENDIX B 

RECOMBINATION EMISSIVITY 
The integration over the H I column density distribution of QSO absorption-line systems 
appearing in the recombination emissivity term (eq. [24]) can be performed analitically over 
ranges of constant rj = iVnell/^Hl and for rj oc l/iVni using approximated formulae for the 
absorption and escape probabilities. Taking 

l-e"^?^ (Bl) 

l + 0.75r ^ ^ 

m the formulae for Pem? ^abs 

and (i, the integral 

fNM Q2j^ 

and the integral over the column density distribution which appears in the expression for the 
cosmic effective optical depth (eq. [5]) can be solved analytically. The computations are simple 
but lengthy. Though the right hand side of eq. (Bl) does not have the right limit for large optical 
depths, it is a better approximation than r/(l + r) for r < 5, with an error for r 1 which is 
lower than 10%. The typical error of the analytical approximations presented in this Appendix 
is of order a few percent, with maximum values < 15%. 



a) Effective Optical Depth The approximation given in eq. (Bl) can be used to calculate the 
integral over the column density distribution which appears in eq. (5). When v < Avl, or when 
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V > 4i^i and the integration is over ranges of costant r\ (i.e. for iVni < iV^j = 1.2 X 10"^^ cm 
we can write 

^(//) = ^iVo(l + ^r^arctanyiVmA , (B3) 

where the constant iVo and the exponent 7 are specified according to eq. (7). We have defined 
r = \jYj{v\ with S(i^) = [3/4r(i^)]/iVHi. When > 4i^i, and the integration is over ranges of 
1] oc l/iVni (i.e., for iVni > ^hi)' have 



1 



1 + A(//) 



arctan iVni / ? 



A(^) 
1 + A(//) 



1 



HI 



(B4) 



where tv{v) = [3/4iVgj?7maxCHell(i^)]- Here r = [1 + A(i^)]/E(i^), with E taking into account only 
the opacity due to H I and He I . 

b) H 1 Recombination Continuum In computing the emissivity due to H I LyC radiation, we 
first note that the contribution of photons with > Avl to the integral over v' appearing in eq. 
(24) is negligible, and that the emissivity at > Avl is small compared to the emissivity due 
to He II LyC radiation at the same frequencies. Thus we can evaluate eq. (B2) for vl < v and 
v' < Avl, to obtain 



2 ^ j^. 

K{v, v') = No(l + zy r2r3(7m(i'') arctan ^/Nm/r^ 
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for V ^ v', where E(z/) = [3/4r(z/)]/iVHi, as = aHell-ai^s, ro = 1/S(1 Ryd), ri = aHell/[S(l Ryd)^^] 
r2 = 1/S(z/'), rs = 1/E(z/), and 



Ri = [ro - ri]/[{r2 - ri){rs - ri)]; 
R2 = [ro - r2]/[iri - r2){rs - r2)]; 
R3 = [ro - r3]/[{ri - rs){r2 - rs)]. 



(B6) 



It is easy to obtain a similar formula when v = v' . In this case the second order multiplicity 
of the root S(i^) leads to a different decomposition of the integrand into integrable elementary 
partial fractions. We find 



K{v', v') =iVo(l + zy — rlumiv') 

2 



— zz arctan Nju/ri -\ ^ -\ ^= arctan iVni / ^2 
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r2Nm + rl r2^Jr2 



Nm (B7) 



where now 



Ri = 


[ro 


- ri]/{r2 
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[ri 


- ro]/ir2 
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- r2]/iri 


- r2). 
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c) H I Lya and Two-Photon Continuum The integral in eq. (B2) can be ulteriorly simplified 
when Pem{v) = 1, e.g., in the case of H I Lya and two-photon emission. We obtain 



Z , \ - JXi 



K{v') = Nq{1 + zy r2(Tm{v') > ^ arctan ViVm/r^ 

1=1 ^ 



where 

Ri = [ro - ri]/{r2 - ri); 
R2 = [ro - r2]/iri - r2). 



(B9) 



(BIO) 



The formulae (B5), (B7), and (B9) above are strictly valid for any integration range only when 
hiy' < 24.6 eV. At higher photon energies, they are only correct as long as Nm < ^Hl(4Ryd), as 
clouds with larger columns are characterized by a larger He I to H I ratio. The correct form of 
WabsiNm^ v') for these thicker clouds is 

WaUr, v ) = const + [1 - e'^^^')] ^^i^. (BH) 

The integral over column in eq. (B2) can then be split into two distinct terms and solved. 

d) He II Recombination Continuum The evaluation of the emissivity due to He II recombination 
continuum can be simplified noting that the absorption of He II ionizing photons occurs at the 
He II photosphere. We also note that photons above 54.4 eV are effectively absorbed only by 
He II . The integral in eq. (B2) for iVni > iVHi(4Ryd) can then be written as 

Ki^, = ^ Pe^(^)%^(^'), (B12) 

a A — U.o I ai2s az 

where pem is given by eq. (28). The result for iVni < iVHi(4Ryd) is analogous to eq. (B5), with 
(Thi replaced by (JHell and E evaluated at 4 Ryd in ro and ri . 

e) He II Lya and Two-Photon Continuum The arguments given above also imply that, in 
the case of He II Lya and two-photon emission, pem — 1- The integral in eq. (B2) for iVni > 
iVHi(4Ryd) can then be written as 

K(iy') = ^^(//'). (B13) 

^ ^ aA- 0.37ai25 dz ^ ^ 

This becomes analogous to eq. (B9) for iVni < iVHi(4Ryd), again with cthi replaced by (JHell and 
E evaluated at 4 Ryd in ro and ri . 

APPENDIX C 

NUMERICAL INTEGRATION OF THE EQUATION OF RADIATIVE TRANSFER 
The key problem to address for an iterative solution of the radiative transfer equation is how 
to achieve a fast numerical estimate of the effective optical depth Tg// in each iteration. This 
is a three-dimensional function, since at any observed frequency it depends on the observer 
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and emission redshifts Zq and z^m- Computing time can be optimized noting that in the exphcit 
form of Tg// (eq. [5]), these three variables appear in the combination v = Vo{^ + ~\~ ^o)- 

A logarithmic redshift grid in the variable x = 1 -\- z is therefore designed, so that the ratio 
R = Xj^xjxj is constant for any j < J. We also set up a frequency grid longward of the H I and 
He II Lyman edges, defined by 

v, = vlZ'^ IR'-^ I = 1,2, ...I. (CI) 

We make sure that the Lya frequency coincides with one of the sampling points in this grid by 
properly adjusting the value of i?, i.e., by setting R = (4/3)"^/"^, where m is an integer > 1. A third 
frequency grid provides coverage of the frequency range above 54.4 eV. The next step consists of 
computing the two sets of integrals which are needed to evaluate Tg// as well as dreff/dz. The 
first of them, 

where 

Tzj = NmWmii'iXz/xi) + r]{xj)aiieiiii'iXz/xi)], (C3) 
is computed for 1 < z < / and j > i. A second set of integrals, 

where 

Tzj = NmWmii'iXi/xi) + r]{xj)aiieiiii'ixi / x,)], (C5) 

is evaluated for 1 < z < / and 1 1^ j 1^ J ■ Note that U(l,j) = L(l,j), and L(i,j) = longward of 
the H I Lyman limit. We can now calculate the effective optical depth re//(i^o, -^o, -^em) = T(i,j, k) 
(i ^ k) over the frequency X redshift X redshift space. Using the fact that the redshift distribution 
of the clouds can be expressed as a power-law (eq. [7]), we find the following recursive relation: 

T(l, J, k) = T{l,j, k-l) + [R^ik -J + l,k) + Uik-j,k- l)]dk-i (k > 2); 

T{i,j,j + 1) = [R-^Lii + + L{i,j)]d, (z > 2); (C6) 

T{i,j, k) = T{i,j,j + 1) + T(z - l,j + 1, k), 

where 

d.=xr^^-^. (C7) 

Note that T(i,j,j) = 0, and T(i,j, k) = for k < i -\- j — 1 longward of the H I Lyman limit. 

A computer code, named CUBA (from Cosmic Ultraviolet BAckground), has been designed 
based on the scheme sketched above. For each of the three frequency ranges, CUBA computes 
a number of integrals less that 2(1 X J). It is worth to point out that one should be careful in 
estimating the cloud emissivity (as described in Appendix B) at the photoionization threshold 
frequencies. Indeed, the opacity for photons at observed frequency < vi is zero for Zg < 
z < {vilvj)zo, while the volume emissivity is non-zero for z > {vilvj)zo. The spectra shown 
in Figure 5 have been computed by CUBA over 35x165 time X frequency grid points in a single 
run. Convergence (1 part in 1000) was reached after 16 iterations. Each iteration took about 15 
seconds of CPU on a Digital VAX Alpha 3000/300. 
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Figure Captions 

Figure 1: Solid curves: The He ll /H l ratio as a function of H I column density through a 
slab illuminated on both sides. The assumed total hydrogen gas density and temperature are 
2 X 10""^ cm~^ and 2.5 x 10"^ K, respectively. The three curves refer, from top to bottom, to a meta- 
galactic flux at the Lyman edge of intensity 5 X 10~^^ , 10~^^ , and 5 X 10~^^ ergs cm~^ Hz~"^ sr~"^ . 
The spectral index is taken to be a = 2. Dashed curves: Our "step-like" approximation (see text 
for details). Dotted curves: same ratio obtained by self-consistently solving the thermal structure 
of the clouds. 



Figure 2: Solid curves: Absorption probability (1 — e~'^) through an individual cloud as a 
function of neutral hydrogen column for four different ionizing photon energies (54.4, 75, 100 and 
300 eV). Dashed curves: Same with our step-like approximation to the to the probability. Dotted 
curves: Same assuming a constant He II /H I value as derived in the optically thin limit. All 
curves are relative to the J912 = 5 X 10~^^ ergs cm~^ s~"^ Hz~"^ sr~"^ background model adopted 
for Figure 1. 

Figure 3: Upper panel: Average temperature of a cloud as a function of neutral hydrogen column 
for the three background models adopted for Figure 1 (solid line: 5 X 10~^^; dotted line: 10~^^; 
and dashed line: 5 X 10~^^ ergscm"^ Hz~"^ sr~"^). The gas total hydrogen density is taken to 
be 2 X 10""^ cm~^ for iVni < 1.59 x lO^^^cm"^, and 4 x 10~^ cm~^ for thicker systems. Lower 
panel: Fractional population density of the n = 2 H I level due to recombinations of photoionized 
ions (upper curves), coUisional excitations of neutral atoms (middle curves) and recombinations 
of coUisionally ionized ions (lower curves), as a function of the neutral hydrogen column. 

Figure 4: The QSO proper volume emissivity at (from top to bottom) the H I , He I and 
He II ionization edges as a function of redshift. See text for details on the assumed QSO lu- 
minosity function and spectral energy distribution. 

Figure 5a: The ionizing background from 5 to 5000 A (solid lines) at epochs z=0, 0.5, 1 and 
1.5, computed in a go =0.1 cosmology taking into account the absorption and reprocessing of 
QSO radiation by discrete absorption systems. The spectrum is compared to that arising from 
a perfectly transparent (te// = 0) universe (dotted lines) ^ and a purely absorbing intervening 
medium (dashed lines). 

Figure 5b: Same as Figure 5a but for z=2, 2.5, 3 and 3.5. 
Figure 5c: Same as Figure 5a but for z=4, 4.5, 5 and 5.5. 
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Figure 6: Upper panel: Photoionization rates for H I (solid line), He I [dashed line) and He II [dot- 
dashed line) derived from our model background as a function of redshift. Lower panel: Ratios 
between these rates and the photoionization rates computed in the hmit of a purely absorbing 
intervening medium. We also plot (again relative to the same quantity in the limit of a purely 
absorbing medium) the ratio between the O II and C III photoionization rates. As the O II and 
C III ionization potentials are 35 and 48 eV, respectively, this quantity will be particularly sen- 
sitive to He II Lya line and two-photon continuum recombination radiation from intergalactic 
absorbers. 



Figure 7: Solid curve: The He ll /H l ratio as a function of redshift for an optically thin gas in 
ionization equilibrium with the ultraviolet background shown in Fig. 5. Dashed curve: Same in 
the limit of a purely absorbing medium. 



